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Supported by NASA Grant NAG-1 -4 16 
Report Period: 1/1/85 - 12/31/85 

The Solid State Physics Research institute at Virginia State University 
was formally orgauiized in January 1984 subsequent to the funding received 
through NASA grant The institute is a c D ll e ct ie on of three- research 

programs. Two of these, muon spin rotation studies and studies of 5mnealing 
problems in galliion arsenide, were funded previously by NASA through separate 
grants, ihe MuSR program was beg un in 1972, while l±e GaAs studies were initiated 
in 1982. A third program. Hall effect studies in semiconductors, was initiated 
with the establishment of the institute. C. E. Stronach is director of the in- 
stitute as well as P/I of the MuSR program. J. J. Stith and J. C. Davenport 
(VSD phTSiea dapactaSBit ebsizan) are co-principal investigators of the GaAs 
program, «^iile G. W. Henderson is P/I of the Hall effect progr a m. 

The purposes of this institute are threefold: (1) to perform state-of-the- 
art research in both basic and applied aspects of solid state physics which is 
both intrinsically interesting and which can be applied to problems of i n te r est 
to NASA; (2) to develop a self-sustaining physics research program at vsn wh ic h 
will bring djatinction to the university and its physics departm^it; and (3) to 
provide training and experiences to students which will prepare them for careers 
in research. 

The following sections describe tlie activities of each of the research 
programs during the period Janueury 1 - December 31, 1985. Following this is a 
section describing activities which do not fall strictly within any of the three 
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programs# plus descriptions of student activities# equi^ent acquisitions# 
and future plans. 


Muon Spin Rotation .Studies 

Ibe ai»n span rotation project efforts centered arosoad experinents done 
at the Alternating Gradient Synchrotron of Brookhaven National Laboratory, 
jageriaents and beam development done in 1983 and 1984 were followed by a 
number of experiments done in March/April emd May /June 1985. These include 
studies of TlH,^ and YH 2 # TMMC# aluminum alloys# and heavy fermion supercon- 
ductors. Appendix 1 is an edastracb of a paper to be presented at the March 
1986 American Physical Society meeting on the TlH^ results# along with a plot 
of the date* Appoaddhr-2 isvem adbetsest Hor en amalogons peper on the aiuninam 
3d.loy data. 

Analysis of iron alloy data shows that the magnetic field sensed by the 
muon# Bp# is greater in magnitude for Fe(Au) than for pure iron# but it is 
snaJler on magni tade for fe (Dy) and Pe (Ta) . The effect is valence axui hence 
electron density related and follows the "Doyama rule" which was first stated 
after analogous positron studies of alloys. That is# {and e*^) are repelled 
by impurities to the right of iron in the periodic table, but are attracted by 
impurity atoms to the left of iron in the periodic tadsle. The Fe(Ta) results 
are strongly dependent on sample preparation cmd a complete report on these 
materials aw 2 d.ts a study of Fe(Ta) alloys to be done at Brookhaven in May 1986 
in which the shifts of By upon annealing are analyzed in detail. An abstract of 
a paper presented on this topic at the March 1985 American Physical Society is 
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included as Appendix 3. 

We are most pleased to report that the paper, "Uniaxial stress Induced 
symmetry breaking for muon sites in Fe" was published in the July 1, 1985 
issue of Physical Review B. This study was a major ongoing effort for about 
five years (1980-85) and was supported both by NAG-1-416 and the preceding 
grant NSG-1342. A reprint of this paper is included in this report as Appendix 4. 

The principal i awastigatur and two graduate stocbents, Incien Goode* Jr. 
and Nana Adu, visited the Saclay laboratory in France in June/July 1985 to 
study muon spin rotation in Ni and Ni(Ft) in conjunction with the French group 
led by Dr. Robesrt 1. Gbrynszpen of the Oeater for the Stady of Metallurgical 
Chemistry, Vitry-sur-Seine , France. The travel funds were provided by the 
USA/France Muon Spin Rotation Program, which is supported by the National 
Science Foundation. Mo NASA funds were used. The purpose of the experiment was 
to detaxttine whe th eg a'*' ia »t<Ft) trap at sites near the Pt atoms for which 
is substeurtially different fro m that for pure Mi. An extensive study at 77K 
showed no such trapping. 

The principal investigator served on the scientific advisory committee of 
the European Morkshop on Subatomic Spe c ies in Solids, which «ias h eld in Vitry, 
France, September 2~-€, 1985. He also gave a talk on "Qu£uitom field theoretical 
aspects of muon spin rotation** at the meeting. Again, travel was supported by 
the NSF grant for the USA/France Muon Spin Rotation Program amd no NASA funds 
were used. A digest of the talk is included here as Appendix 5. 

Another round of ySR experiments is planned for May/ June 1986. These are 
expected to include additional Fe(Ta) measurements, additional heavy fermion 
superconductor studies, fatigue studies in metals (in collaboration with our 
French colleagues) , and studies of additional metal hydrides. 
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An experiment recently done at the KEK laboratory in Japem (and even more 
recently improved upon at the TRIUMF laboratory in British Columbia) bodes well 
for future ySR research ^ particularly with regeurd to its utilization in the 
study of surfaces. A low-energy y'*' beam was projected onto tungsten foils. Three 
to five percent (20 to 25% at TRIUMF) of the incident y'*' were subsequently de- 
sorbed with enezgies of a few electron volts. Me plan to design a beam based 
this phenomenon to provide extremely low energy y"** at the Alternating Gradient 
Synchrotron (Brooklmcven} for studies of surfaces. Preliminary studies wall 
ccanmence during the May /June 1986 run. 

Prepara t ions tor the 198£ raa JjacJuds deioalopnent of the PDF 11/73 costputer 
(programming, purchasing peripherals and softwaure) and a cryogenics system for 
low-temperatxire studies. Electronic timing should be improved with the acquisi- 
tion of a higher resoltction tine-to-amplitude converter and a co n s t ant-fraction 

The Virginia State University participants in these experiments were Carey 
£. Stronach, director of the institute* Lucian R. Goode* JT. , and Nana Adu, 
physics graduate students, and Michael R. Davis, an undergraduate physics major. 
The Brookhaven pr ogr am is under the overall leadership of William J. Kossler* 
pr ofess or of physics at the College of William and Mary. Other participants 
included William F. Lankford of George Mason University, Harlan Schone of 
William and Maury, Gabriel Aeppli of Bell Laboratories* Alain Taooanc of CNR5, 
Grenoble, France, Ola Hartmann and Roger WSppling of the University of Uppsala, 
Sweden, Joseph Budnick of the University of Connecticut, and two graduate 
students from HiUiam and Mary. 


Computer Simuiat'ion of Radiation Damage 

ORIGINAE PAGE IS 
OE POOR QUALITY 

During the past year three graduate students and one 
undergraduate participated in the radiation damage research 
project. A continuation of this level of student participation is 
anticipated for the coming year. John Stith, one of the principal 
investigators, attended two professional meetings during the 
year<The 1985 IEEE Annual Conference on Nuclear and Space 
Radiation Effects in Monterey, CA and The Eighteenth IEEE 
Photovoltaic Specialists Conference in Las Vegas, Nevada). A 
graduate student, Larry Brown, accompanied John Stith to the 
Nuclear and Space Radiation Effect Conference. A Late News Paper 
which was co-authored by John Stith of Virginia State University 
and John Wilson of NASA Langley Research Center was presented at 
the Photovoltaic Specialists Conference by John .Stith. John Stith 
also spent.^ approKimately six weeks at the NASA Langley Research 
Center during the summer. 

Graduate stude^ participation is an essential component of 
the research grants presently, three graduate students are 
involved in research projects directly related to the radiation 
damage studies. These students will eventually use their findings 
or results to write theses which they will use as a part of their 
requirements for the Master of Science degree. One of these 
students is scheduled to graduate in May. His research project or 
thesis topic deals with an important aspect of radiation damage 
studies utilizing computer simulation studies. He has made a 
comparative study of the effect of using two different potential 
functions in the model used for the simulations. This study is an 
extremely significant aspect of radiation damage studies 
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utilizing computer simulation techniques since the interatomic 
potential is one of the most important parameters to be chosen in 
developing simulations models. 

This aspect of the project demanded a great deal of time and 
input from the principal i nvest i gators who are associated with 
the radiation damage studies. Formal courses on radiation damage 

studies are not offered in the department, consequently the 

/ 

principal investigators must provide individualized instructions 
and supervision for these students participating on the project. 
Consequently, the project has served an important and informal 

t 

mechanism for enhancing students' graduate training and 
experiences. 

The radiation damage simulation research has been 
concentrating mostly on simulations of radiation damage produced 
in a gallium arsenide crystal when it is irradiated with protons 
and electrons. The simulations were performed with a modified 
version of the binary collision simulation code called MARLOWE. 

The modified version -wl||^ch has been described in an earlier 

»*• 

report is referred to as MARsJ Approximately one hundred 
different computer runs were made during the year in order to 
build up a substantial data base and to test the substantial 
number of modif ications that had to be made in the code to keep 
it compatible with the computer operating systems at NASA. These 
modifications were made necessary by the seemingly continuous 
upgrading of the computer system at NASA. In order to perform 
useful simulations with this code an acceptable and realistic 
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model was first developed for the gallium arsenide crystal. One 
of the most important aspects of the model is the selection of an 
interatomic potential function to be used for the interatomic 
interactions that are simulated. The Moliere and the Born-Mayer 
potential functions were considered because of the behavior of 
these two functions for relatively large nuclear separati on (=>1 
A*-’) where the outer electron shells become important. These large 
nuclear separations are expected during collisions between a . 
relatively low energy primary recoil atom and other atoms inside 
the crystal . The Moliere potential function is ' 

V(r)=( (Zi*Z2*e2 )/r)*f(r) 

f (r) = (0.35*exp (-0.30*r/a> +0.55*exp (-1 . 2*r /a) +0. 10*exp <-6.0*r/a) ) 
where 'a' is the screening radius and 'r' is the interatomic 
separation. The Born-Mayer potential function is 
V(r) = Aexp(-rZag) where 'r' is 

the interatomic separation, 'ag' is the screening radius -and the 
average value of 'A' for gallium arsenide is 9,413 eV. Plots of 
the Moliere and Born|^Mayer potential functions for a range of 
interatomic separations are shown in figure 1. The Born-Mayer 
potential function was finally selected as the function to be 
used in the model because of the way in which it approaches zero 
at relatively large interatomic separation and the fact that it 
was derived from experimental results. 

The simulation code has been used to generate results for 
radiation produced damage in gallium arsenide by keeping track of 
the action of a simulated primary recoil atom as it moves among 
the lattice sites of the crystal. The primary recoil atoms are 
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given a range of energies correspondi ng to some of the possible 
energy transfers from 1-MeV electrons to atoms of gallium or 
arsenic. The maximum transfer from the 1-lieV electron is 
approximately 65 eV. The 1-MeV electron was chosen because of its 
importance in equivalent electron fluence studies^. However , 
further studies performed during the year on the concept of 
equivalent electron fluence yielded results that indicated that 
1-MeV electrons should be replaced by 10-MeV electrons in order 
to make the equi val ent ■ el ectron fluence concept more valid. The 
studies are described in greater detail in reference' 4. A copy of 
reference 4 is included in the appendix of this report. 

The results from the simulations include the followings a. 
number of displaced atoms; b. number of improper replacements; c. 
number of improper replacements; d. number of proper 
replacements; e. inelastic energy loss; f. range of the primary 
recoil atom; g. number of interstitials atoms; h. number ’ 
vacancies; and i. number of Frenkel pairs. 

A theorectical Understanding of displacement effects in 
semiconductors requires a detailed analysis of the dependence of 
the quantities listed above on the energies of the primary recoil 
atoms. This analysis is done by plotting graphs and studying the 
relationship between the various results from the simulated 
radiation and the energy of the primary recoil atoms. These 
graphs represent some interesting and important results. A graph 
of the total number of displacements versus energy in figure 2 
indicates an approximately linear relationship between the total 
number of displacements and the energy of the primary recoil atom 

i 
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for the range of energies used in the simulations. Curve 2 in 
this figure was plotted using results obtained from a 
computational model developed earlier by Wilson, Stith, and 
Stock'*. Curve 3, included in this figure, represents a plot of 
the number of Frenkel pairs versus energy. Figure 3 is a plot of 
the inelastic energy losses of the primary recoil atoms versus 
the initial energies of the primary recoil atoms. Figure 4 is a 
plot of the ratios between the inelastic energy losses and the 
initial energies of the primary recoil atoms versus the initial 
energies of the primary recoil atoms. From the graph* in figure 4 
it appears that the percentage of the primary recoil atoms' 
energies that is lost during inelastic collisions approaches a 
maximum of approximately 20 percent. Figure 5 is a plot of the 
total number of defects versus energy. The total number of 
defects is the sum of the number of improper replacements, the 
number of interstitial atoms, and the number of vacancies. 

Further studies will be made on the above results during the 
coming year. A stud>^of the effect of temperature on the 

:V 

production of defects will be included. 
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During the year the Hall effect research prograM woved 
through several nilestones. The equipment was acquired which 
will allow conventional Hall effect measurement to be made at 
yirginia State University. The program involved two 
undergraduates and one graduate students in the system 

development. A working relationship was developed with the Naval 
Research Laboratory was also commenced. This contact will allow 
the performance of Quantum Hall effect. 

We also consulted Or. Robert Coleman at the University of 
Uirginia. He is during Hall Effect experiments. The principal 

investigator also attended a short course in Radiation Damage 

through the auspices of the IEEE a Monterey, California. 

The measurement system at Virginia State University ' 

f 

consists of a Walker Scientific Power supply Cfrom 0 to 50 amp) , 

an alpha Magnet, a digital Gaussmeter, a HP 98tB computer with 

peripherals, a 3497A data acquisition will be used to acquire and 
analyze the data. A cryostat has also been acquired. 

Roscoe Ledbetter, an undergraduate student has 
developed a computer program for data acquisition. A program to 
process and analyze the data is nearing completion. The decision 
has been made to operate two systems of sample geometries. One 
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employs the standard Hall sample configuration. The other will 
utilize the Van der Pauw technique. 

The experimental procedure as we have perceived it 
consists of the measurement of the Hall coefficient, resistivity, 
and Hall mobility as a function of temperature of a sample of 
Gallium Arsenide before and after irradiation with low and high 
energy protons. In order to compare the results with the known 
damages that are produced by proto irradiation , one need to know 
the relationship among the quantities that radiation damage and 
the Hall effect measurements have in common. 

It is known that the effect of radiation damage is 
associated with the carrier concentration and the Hall 
mobilities . The manifestation of the radiation defects 

I 

are "intrinsic defects". In order to establish and test the 
theory , a theoretical model is being developed using 
Lindhard’s theory of atomic collisions. 

The second most important step is to link 
the quantum Hall effect and the changes 
produced by proton collisions. Work has 
begun on the the development of a theory and 
subsequently experiments will be designed and 


per f ormed . 



APPENDIX I 

RESISTIVITY AND HALL COEFFICENT MEASUREMENTS 
VAN DER PAUW METHOD 


The neasureMent theory deve^ioped by Van der Pauw involves the use of 
arbitrary geonetric configurations . Incunbent adjustments are made for the 
chosen geometry. For our system consider a lamella 

In order to make resistivity measurements a current I is passed through two 
adjacent contacts e.g. <2 and 3), the voltage V is then measured across the 
other contacts. The equation R == V/I is used. Next a current I’ is passed 
through the next pair of contacts (3,4) and the voltage V’ is measured across 
the other contacts. The equation R’ « VVI’ is used. The resistivity * ’ 
of sample of thickness ’t’ is related to R and R* through the relation 
Exp(- tR/ ) + Exp<“‘ tRV ) * I, The resistivity can thus be obtained. 

In order to determine the value of the Hall coefficient, a current I is 
passed through two non adjacent contacts and the voltage V is measured across 
the other two. The relation R ® R/I is used. The magnetic field is now 
energized and the mea|prements of V* and I’ are then made on the same 
contacts. The Hall coefficient can now be determined from u « R^./ 

The measurement process consists passing a current through both 

directions and averaging the results. The magnetic field is reversed and 
measurements taking in both direction and averaged. 

The following relations are typical for the experimental determination of 
the resistivity and Hall mobility. All measurement are made as a function of 
temperature . 

- t/ln2(V„,/2I,, ^ ) f 

23 14 24 31 

Uu - )21n2/ 

H 32 21 34 13 42 



10 
1 1 


CURVE FITTING PROGRAM FOR HALL EFFECT DATA HENDERSON JORDAN FEB 7 . 1 98G 
IDENTIFIERS 


13 

! Cl 

COEF 

SOLUTION VECTOR 

14 

! C3 

CORREL 

CORRELATION COEFFICIENTS 

15 

! E2 

SIGMA 

VECTOR OF ERRORS 

16 

1 E5 

SEE 

STD ERROR OF ESTIMATE 

17 

! L3 

NLIN 

NUMBER OF PLOT LINES 

18 

! Ml 

MAX 

MAXIMUM LENGTH 

19 

! Nl 

NROW 

NUMBER OF ROWS 

20 

i N2 

NCOL 

NUMBER OF COLUMNS 

21 

! R5 

RE5ID 

VECTOR OF RESIDUALS 

0 7 

! 57 

SUMY 

SUM OF Y 

23 

• SB 

SUMY 2 

SUM OF Y SQUARED 

24 

1 T6 

SRS 

SUM OF RESIDUALS SQUARED 

25 

! Y3 

YCAL 

CALCULATED Y 

26 

i END OF IDENTIFIERS 


40 

A$= '■ DDO OD . 0 " 


50 

C*="DDD DD.d" 


60 

f^l=35 



70 

DIM Z<4).A{4.4),C1< 14). Y( 

35>.U<35,4) ^ 

80 

DIM W<4,1 ).B<4, 4). 12(24,2 

3>,X(35).Y1\35) C) 

90 

DIM Y2(35).R3<35),E2<4) 


100 

GRAPHICS OFF 


1 10 

PRINT 



120 

PRINT “ UP 

TO THIRD DEGREE POLYNOMIAL LEAST SQUARE’ 

LIMINATION" 



130 

PRINT 



140 

605UB 500 

! GET THE 

DATA 

150 

! SORT THE 

DATA 


160 

GOSUB 800 

! SQUARE 

UP THE MATRIX 

170 

60SUB 4000 

! SET UP 

THE MATRIX 

180 

GOSUB 5000 

! GAUSS- 

JORDAN SOLUTION 

190 

GOSUB 1000 

1 PRINT 

THE RESULTS 

200 

GOSUB 7000 

! PljiT DATA 

210 

GOTO 100 



500 

! GET THE 

DATA 



FIT BY GAUSS JORDAN E 


510 

520 

530 

540 

550 

560 

570 

580 

600 

610 

620 

630 

640 

650 

660 

800 

810 

820 

830 

840 

850 

860 

870 


INPUT "NUMBER OF DATA POINTS". Nl 
INPUT "POLYNOMIAL ORDER" .N2 
IF N2>3 THEN 520 
IF N2C1 THEN 9999 
N2=N2+1 
L3=(N1-I )*2+l 

IF N4=1 THEN 620 
FOR 1=1 TO Nl 
READ X<I>,Y1(I ) 

NEXT I 
RETURN 
! Y DATA 

DATA 1 .2.07,2,8.6.3,14.42,4.15.8,5,18.92.6 .17.96 
DATA 7, 12. 98. 8. 6. 45. 9. .27 
RETURN ! FROM INPUT ROUTINE 
! SET UP THE MATRIX DATA 
FOR 1=1 TO Nl 


U(I .1 )=1 
FOR J=2 TO 
U< I ,J >=U< I 
NEXT J 
Y< I )=YH I ) 
NEXT I 


N2 

J-1 )*X( I ) 



DOW rvciuniN ! rnuii isc.ii.ii'to v< uiim iduima 
1000! CALCULATE RESIDUALS AND PRINT RESULTS 
1001 PRINTER IS 701 
1010 S7=0 

1020 S8=0 

1030 TG“0 
1040 FOR 1=1 TO N1 
1050 Y3=0 

1060 FOR J=1 TO N2 
1070 Y3=Y3+C1( J )*U( I ,J ) 

10B0 NEXT J 
1090 R3(I)=Y3-Y(I) 

1100 Y2(I)=Y3 

1110 T6=TG+R3< I)''2 

1120 S7=»S7+Y<I) 

1130 S8=S8+Y<I)"2 

1 140 NEXT I 

1 150 C3=SQR< l-T6/<S8-S7'-2/N1 )) 

1160 IF N1=Ns THEN E5=SQR< T8 ) 

1 170 IF N10N2 THEN £5=SQR< T6/< N1 -N2 ) ) 

1 1 B0 FOR J = 1 TO Nil 
1 1 90 E2U )=E5*SQR< ABSl B< J .J))) 

1200 NEXT J 

1210 PRINT '■ X Y YCAL RESID" , 

122G FOR 1=1 TO N1 

1230 PRINT USING "5X.DD,6X, DDD.DD .5X ,0DD = D0 ,6X .ODD.DD .6X .DDD.DD" ; I ,X< I ) ,Y< I ) ,Y 
2(I),R3<I) 

1240 NEXT I 
1250 PRINT 

12B0 PRINT “COEFFICIENT ERRORS" 

1270 PRINT USING ''5X JJDD.DD .3X ,DDD.DD'‘ ;C1< 1 ) ,E2( 1 ) 

1230 FOR 1=2 TO N2 

1300 PRINT USING “5X .DDD.DD .3X ,DDD.DD” 5 CI < I ) ,E2C I ) 

1310 NEXT I 
1320 PRINT 

1330 PRINT USING " 1 0X ,27A . 1 X ,DDD . DO” s “CORRELAT ION COEFFICIENT IS” .03 

1331 PRINT 

1332 PRINT 

1333 PRINT "EQUATION OF CURVE " 

1334 PRINT USING " 1 0X .ODD . DO . 1 X , 1 A . 1 X ,DDD. DD . 1 A , 1 X . 1 A , 1 X .ODD . DD .3A . 1 X , 1 A .ODD .00 

,5A" iCI <1 ) " ,01 (2 ) ."X" .‘'+#,01 (3 ) ,"X*X” ” .C1< 4 ) , "X*X-*X " 

133G WAIT 10 

1337 PRINTER IS 1 

1340 RETURN 1 FROM PRINTOUT 

4000 ! U AND Y CONVERTED TO A AND Z 

4011 ! IDENTIFIERS 

4013 1 N1 NROW NUMBER OF ROWS 

4014 ! N2 NCOL NUMBER OF COLUMNS 

4015 ! END OF IDENTIFIERS 

4020 FOR K=1 TO N2 

4030 FOR L=1 TO K 

4040 A<K,L)=0 

4050 FOR 1=1 TO N1 

4060 A<K ,L>=A<K .L )-t-U< I ,L )*U< 1 ,K > 

4070 IF K<.'‘L THEN A'i L ,K >=A( K ,L) 

4080 NEXT I 

4090 NEXT L 

4100 Z< K )=0 

4110 FOR 1=1 TO N1 

4120 Z<K )=Z<K )+Y( I )*U< I ,K ) 

4130 NEXT I 

4140 NEXT K 

5000 ! JORDAN MATRIX INVERSION ANO SOLUTION 

5010 1 FEB 7, 1986 

5011 I IDENTIFIERS 

5012 ! A A COEFFICIENT MATRIX 


ORJGWai; prcr „ 
Quautx 


S.JV t V 


5014 

5015 

5016 
501 7 
5016 

5019 

5020 

5021 

5022 

5023 

5024 

5025 

5026 

5027 
5029 
5060 
5090 
5100 
5110 
5120 
5130 
5 ! 40 
5150 
5150 
5170 
5180 
5190 
5200 
5210 
5220 
5230 
5240 
5250 
5260 
5270 
5280 
5290 
5300 
5310 
5320 
5330 
5340 
5350 
5360 
5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
5480 
5490 
5500 
5510 
5520 
5530 
5540 
5550 
5560 
5570 


B1 

BIG 

BIGGEST VALUE 

Cl 

COEF 

SOLUTION VECTOR 

0^ 

DETERM 

DETERMINANT 

E 1 

ERME3 

1-“ rN r* 1 a r* 

r-H.nun runo 

HI 

HOLD 

WORK VARIABLE 

12 

INDEX 

WORK MATRIX 

13 

IROW 

ROW INDEX 

14 

ICOL 

COLUMN INDEX 

15 

INVRS 

PRINT- INVERSE FLAG 

N2 

NCOL 

NUMBER OF COLUMNS 

N3 

NVEC 

NUMBER OF CONSTANT VECTORS 

PI 

PIVOT 

PIVOT INDEX 

W 

W 

SOLUTION MATRIX 

Z 

Z 

CONSTANT VECTOR 


i END OF IDENTIFIERS 
! 

El-0 ! BECOMES 1 FOR SINGULAR MATRIX 
15-1 ! PRINT INVERSE MATRIX IF ZERO 
N3-1 ! NUJMBER OF CONSTANT SECTORS 
FOR I-l TO N2 
FOR J-! TO N2 
B<IJ;<-A<I.J) 

NEXT J 
W< I ,1 )-Z< I ) 

I2( I ,3)=0 
NEXT I 
03=1 

FOR 1=1 TO N2 
! 

! SEARCH FOR LAR6EST < PIVOT) ELEMENT 


B1*0 

FOR J=1 TO N2 

IF I2<J.3)=1 THEN 5350 

FOR K=1 TO N2 

IF I2<K.3)>1 THEN 6120 

IF I2(K .3)=1 THEN 5340 

IF B1>=ABS<B< J ,K )) THEN 5340 

I3=J 

I4=K • t 

B1=ABS<B< J ,K ) ) 

NEXT K 
NEXT J 
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I2< 14 .3 )-I2< 14 ,3)+1 

I2<I.t)=I3 

I2< I .2 )-I4 

! INTERCHANGE ROWS TO PUT PIVOT ON DIAGONAL 

IF 13=14 THEN 5540 

03=-D3 


FOR L=1 TO N2 


H1=B( 13 ,L) 

6( 13. L )-B( 14. L) 
B(I4.L)-H1 
NEXT L 


IF N3<1 THEN 5540 
FOR L=I TO N5 
H1=W( 13 .L ) 

W{I3.L)-W( 14. L) 

W( 14,1. )=H1 
NEXT L 

! DIVIDE PIVOT ROW BY ELEMENT 
P1=B( 14.14) 

D3-D3-»P1 
B<I4.I4) = 1 
FOR L=1 TO H2 



5590 
’5600 
5610 
5620 
5630 
5B40 
5650 
5660 
5670 
5680 
5690 
5700 
5710 
5720 
5730 
5740 
5750 
5760 
5770 
57B0 
5790 
5800 
5810 
5820 
5830 
5840 
5850 
5860 
5870 
5880 
5890 
5900 
5910 
5920 
5950 
5940 
5950 
5960 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6100 
6110 
6120 
6130 
6 1 40 

7000 

7001 

7002 

7003 

7009 

7010 
7020 
7030 
7040 


NEXT L 

IF'N3<1 THEN 5660 
FOR L=1 TO N3 
W( 14 ,L )=W< 14 ,L )/F) 

NEXT L 

I 

! REDUCE NONPIUOT ROWS 
FOR U = ] TO N2 
IF LI =14 THEN 5770 
T=B<L1 .14) 

B(L1 ,I4)=0 

FOR L=1 TO N2 

B(L1 .L )=B(L1 .L )-B< 14 ,L )*T 

NEXT L 

IF N3<1 THEN 5770 

FOR L»1 TO N3 

W(L1 ,1. )=W(U .L )-W< 14 ,L )*T 

NEXT L 

NEXT LI 

NEXT I 
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j 

FOR 1=1 TO N2 
L=N2-I-H 

IF I2(L,I ) = i2(L.2) THEN 5920 
I3«I2(L,1 ) 

I4=I2(L .2 ) 

FOR K=! TO N2 
H1=B<K .13 ) 

B(K.I5)=B<K.I4) 

B<K ,14 )=H1 
NEXT K 
NEXT I 

FOR K*1 TO N2 

IF 12<K .3)0] THEN 6120 

NEXT K 

E1=0 

FOR 1=1 TO N2 
CH I )=W( 1.1) 

NEXT I 
IF 15= 

PRINT 
PRINT 

FOR 1=1 TO N2 
FOR J=1 TO N2 
PRINT USINS ■' 

NEXT J 
PRINT 
NEXT I 
PRINT 
PRINT ' 

RETURN 
El = 1 
PRINT ■ 

RETURN 




I 

THEN 6140 
MATRIX INVERSE" 


0X .DDD.D" !B< I ,J ) 


DETERMINANT =",03 
! IF INVERSE IS PRINTED 


' ERROR-MATR IX SI NGULAR “ 

! FROM 6AUS5-J0RDAN SUBROUTINE 
PRINT "DO YOU WANT A GRAPH OF THE FUNCTION (Y/N)" 
INPUT K$ 

IF K$=“Y" THEN 7009 
IF K$="N" THEN 9999 

! PLOT OF Y AND Y2 AS A FUNCTION OF X. FEB. 5.1986 
C$=CHR$<255)8r"K" 

DUMP DEVICE IS 701 .EXPANDED 

OUTPUT 2 USING "#,K"tC$ ! Clear leftover display 

PRINT 



7060 

7070 

7080 

7081 
7090 
7100 
7110 
7120 
7130 
7140 
7150 
7160 
7170 
7160 
7 1 90 
7200 
7210 
7220 
7230 
7240 

7250 
72B0 
7270 
7280 
7290 
7300 
73 1 0 
7320 
7330 
7340 
7350 
7360 
7370 
s 

7380 

7390 

7400 

7410 

7420 

7430 

7440 

7450 

7460 

7470 

7480 

7490 

749) 

7492 

7493 
7500 
7510 
7520 
7530 
7540 
7550 
75B0 
7570 
9939 


OUTPUT 2 USING 

PRINT "WHAT X AND Y LABELS DO 
INPUT T1$.T$ 

INPUT "[’lAXlNUM X.Y",N4,N5 
6INIT 

PLOTTER IS CRT ."INTERNAL" 
6RAPHICS ON 

X_odu_max=!00*MAX< 1 .RATIO) 
Y_Qdu_nax=100*MAX( 1 .1 /RATIO) 
L0R6 6 
DEG 

LDIR 90 
CSIZE 4.5 

HOUE 0.Y_Qdu_nflA/2 

LABEL T$ ! Write Y-axi; 

L0R6 4 


! Clear screen for graph 
YOU WANT ON GRAPH" 


Initialize venous graphics parai^etersGRAP 
Use the internal screenGRAPHICS OFF 
Turn on the graphics screenGRAPHICS OFF 
Determine how many GDUs wide the screen is 
Determine how many GDUs 
Reference point; center 
Angular mode is degrees lused in LOIR) 
Specify vertical labels 
Specify smaller characters 
Move to center of left edge of screen 
label 

1 Reference point: center of bottom of label 


high the screen 
of top of label 


IS 


LOIR 0 ! Horizontal labels again 

MOVE X„Qriu_max/2 , . 07*Y_gdu,._,max ! K: center of screen) Y: above key labels 

LABEL 11$ 1 Write X-axis labelGRAPHICS OFF 

VIEWPORT . ! *X_Qdu,_mdx . . 98«X_gdu_max , , 1 5*Y_gdu_max , . 9*Y_gdu_ms(x 

! Define subset of screen area 

FRAME . ! Draw a box around defined subset 


WINDOW 0,N4,0,N5 

AXES 1,1,0.0.10,10,5 ! Draw 

CLIP OFF 

CSIZE 1 .6 ,1 -S 

LQRG 6 

FOR I“0 TO N4 STEP N4/10 
MOVE I 1 

LABEL USING "#,DD.0”iI 
NEXT I 
LORG 8 

FOR 1=0 TO N5 STEP N5/I0 
MOVE -.05,1 


axes with appropriate ticks 
i So labels can be outside VIEWPORT limits’ 

! Smaller chars for axis labelling 
! Ref. pt: Top center 1\ 

! Every 10 units 1 \ 

I A smidgeon below X-axis I > Label X-axis 
! Compact) no CR/LF 1 / 

! et sequens 1/ 

1 Ref. pt: Right center 1\ 

! Every quarter I \ 

! Smidgeon left of Y-axis i > Label Y-axi 


LABEL USING "lt.DD.D"iI 1 
NEXT I 
PENUP 
1=0 

FOR 1=1 TO N I t 

PLOT X( I ) .Y< I ) 

NEXT I 
MOVE 0.0 
FOR K=1 TO N1 
PLOT X(K ) ,Y2<K ) 

NEXT K 

PRINT "IF YOU WANT A HARD COPY OF GRAPH PRESS DUMP GRAPHICS" 
PRINT "IF NOT PRESS CONTINUE" 

PAUSE 

GRAPHICS OFF 

PRINT "DO YOU WANT TO PLOT ANOTHER SET OF DATA Y/N" 

INPUT Z$ 

LABEL TIME$( TIMEOATE ) ,DATE$< TIMEDATE ) 

IF Z$="Y" THEN 100 
IF Z$="N" THEN 9999 
RETURN 

GRAPHICS OFF 
OUTPUT 2 USING "#,K")C$ 

END 

ORIGINAL PAGE IS 

OF POOR quality: 


DD.D) no CR/LF 
I et sequens 


1/ 
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student Participation 


During the Spring 1985 semester six graduate students were participating 
in the program, Lucian R. Goode, Jr., Akpan E. Akpan, Nana Adu, William Bolden, 
Larry Brown amd Karamali Shojaei. Goode, Adu and Shojaei are working on MuSR 
projects; Akpem, Bolden and Brown are working on computer modeling of radiation 
damage in solids. In September 1965 another graduate student, Li-Tai Song, 
joined the Institute and is working on Hall effect studies. 

' Daring the Spr ing 1965 semester f o u r tmdergraduete students participated, 
Michael Davis, Roscoe Ledbetter, Cornelia Belsches and Raymond Noel. Fall 1985 
undergraduates uiare Xkauis, JLedbetter, ReJ aches and Toz^y £arnes. 

Other Activities 

The director of the institute participated in an experiment on inelastic 
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scattering of polarized protons from C at the Los Alamos Meson Physics Facility 
in New Mexico during July /August and November/December 1985. Karamali Shojaei, 
a VSU graduate stxident, participated in the July /August run. Other collaborators 
included Bernard J. Lieb of George Mason University, Herbert 0. Funsten, Charles 
F. Perdrisat and J. Michael Finn of the College of William and Mary, Hans S. 

Plendl of Flora da State university, Jos^h Comfoart of Arizona State University, 
and one greuiuate student each from William and Maary, Florida State ^md Arizona 
State. 

An abstract of a paper to be presented at the April 1986 APS meeting is in- 
cluded as Appendix 6. We are seeking other sources of funding for our partici- 
pation in these eacperiments and hope for success in the near future. 

Two papers based on work done at the Tri-University Meson Facility (Vancouver, 
BC) in 1979-80 and supported in part by NASA grant NSG 1646 were completed in 


1985 with extremely minor support from NAG-1-416. Both were accepted for 
publication by Physical Review C. "Energy dependence of the 'Li(p,d)'’Li 
reaction" appeaured in the September 1985 issue and a reprint is included 
here as Appendix 7. "^He(p,d)^He reaction at 200 and 400 MeV" is scheduled 
to appear in the February 1986 issue and a preprint is included here as 
Appendix B. 

The director participated in the revision of a paper baaed on pion?- 
nucleus reaction experisients conducted at LAMFF in 1960 and supported in 
part by NASA grant NSG-1646. This revised version has been submitted to 
Physiocd. Review c axad a prcjarini- ia j nrJ »ded here as Append j x 9. 

The director of the institute has served as the Virginia State Univer- 
sity representative on the boeurd of trustees of the Southeastern Universities 
Research Association since October 1983. During 1985 he served on the SURA 
industrial affiliates ccnmittee and in early Jaxuiary 1966 was appointed to 
the science amd technology committee of SURA. This committee will provide 
SURA oversight for the Continuous Electron Beam Accelerator Facility, vdiich 
is a 4-GeV continuous wave electron accelerator to be constructed in Newport 
News, Virginia. 

In Hay 1985 the director completed a one-year term as chairman of the 
a sti o no Biy, physics section of the Virginia Academy of Science. 

He was also appointed to the local organizing committee for the International 
Symposivim on the Physics and Chemistry of Small Clusters, which will take 
place in October 1986 in Richmond, Virginia. 

James C. Davenport served as director of the sumner student program at 
Fermilab (Batavia, Illinois) during the summer of 1985. He also served on the 
Committee on Minorities in Physics of the American Physical Society this past 


year . 


George W. Henderson, John J. Stith euid Larry D. Brovm (graduate student) 
received course credit in the summer of 1985 for an lEEE-NSRE tutorial short 
course on radiation effects through the New Jersey Institute of Technology. 


Equipment and Supplies 

The following items were pturchased during the reporting period: 

1 Janis supertran helium transfer tube 
4 EMI 9907B photomultiplier tubes 
3 Ortec 265 tube bases 
3 Ortec 218 phototube shields 
1 EG&G/Ortec 567 Time-to-an^litude converter 

1 PDF JJL/73 co apuf g (l^megabyte memory, 31-megabyte bard diaJc drive, 
dual floppy disk drives, terminal printer 

1 Interface Stemdards IS-ll/CC CAMAC crate controller 

2 Keithley 175 multimeters 
1 Keithley 197 multimeter 

3 Hewlett-Packard HP-15C calculators 

100 reprints of our paper on MuSP in strained single crystals of iron 
1 1-megabyte memory board for HP 9816S computer 
1 Hewlett-Packard terminal emulator 
1 Hewlett-Packard 9816S computer with accessories 
1 Hewlett-Paclrard HP 3497A data acquisition control xaiit 
1 MICRO/RSX operating system for PDP-11/73 computer 
3 PVC inserts and magnetic shields for phototubes 

In addition, the internal account at Brookhaven National laboratory sms 
continued in force. This enedsles the institute to purchase equipment, supplies 
and materials from fiNL directly while experiments are in progress. 


Three peripherals for the POP 11/73 computer were ordered (plotter, modem, 
line filtmH received in early 1986. Several others (IEEE bus, yNlVEWffilt, 
and Advanced Programmers Kit) cure still on order. A quadruple constant fraction 
discriminator and a scientific word processing program are also on order. 



The second year of support from NASA for the Solid State Physics Research 
Institute was a year of growth amd consolidation. The MuSR program saw a major 
publication, substantial progress on several other projects and near completion 
of the new data acquisition etnd analysis system (which should be complete in 
time for experiments in the spring of 1966. The radiation damage studies have 
gone into production mode, and the Hall effect appauratus is essentially ccmt- 
plete. The number of student participants increased substantially over the 
1984 level. 

The higher level of student interest has enahl«?d the principal investiga^ 
tors to set higher stemdards for student eligibility for research stipends, 
and the quality of student involvement showed marked improvement as we begem 
the spring 1986 semester. 

we look forward to another active and successful year in 1966, and we 
appreciate the sup p or t from NASA which makes these activities possible. 


Respectfully submitted. 





Carey E. Stxonach 
Director 

Febrtiary 14, 1986 
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ABSTRACT 


Nuclear 7 -rays in coincidence with outgoing pions or protons 
following single nucleon removal from by 200 MeV have 

been detected with Ge(Li) detectors. Differential cross sections 
are reported for 7 -rays from the first excited mirror states of 
2 ^Na and 23ng in coincidence with positive pions or protons 
detected in particle telescopes at 30*, 60*, 90*, 120* and 150*; 
angle- integrated absolute cross sections and cross section ratios 
<7(23Mg/a(23Na) are calculated. These results are compared with 
the predictions of a Pauli-blocked plane-wave impulse approxima- 
tion (PWIA) and the intranuclear cascade (INC) and nucleon 
charge exchange (NCX) reaction models. The PWIA and the INC 
calculations generally agree with the angular dependence of the 
experimental results but not the absolute magnitude. The NCX 
calculation does not reproduce the observed cross section charge 
ratios. 


PACS number; 25.80 HP 



I . INTRODUCTION 


The (ff,jrN) nuclear reaction offers the potential of a new 
means of understanding nuclear structure and, at the same tiine, 
of studying the propagation of a strong baryon resonance, the a, 
in the nuclear medium. Before this potential can be realized, 
however, a better understanding of the reaction mechanism must be 
developed. 

The ( 7 T,ffN) reaction has been extensively studied using two 
general types of experiment. In one type of experiment, the 
residual nucleus or specific states of the residual nucleus are 
identified through radiochemical techniques [e.g. Ref. 1,2], or 
via detection of prompt de-excitation 7 -rays [e.g. Ref. 3]. 

Since no kinematical or angular information is obtained, these 
experiments integrate over both quasifree and non-quasifree 
components . 

Interest in these experiments has centered on the charge 
dependence of the cross section ratios which were found to differ 
from the free wN ratios. For example, the cross section ratio 
<r(ir'^)/a{n“) for n induced neutron removal from was found to 
have a value of 1/1.7 instead of 1/3 [1]. 

In the other type of experiment [e.g. Ref. 4,5], the 
outgoing pion or proton is detected. With the proper geometry, 
the experimenter can select quasifree events; but only recently 
has the charged-particle energy resolution become sufficient to 
identify the final nuclear state. In the work of Kyle et 
al . [5], the ir'^/ir~ ratio for proton removal from was found 
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to be as large as -40 at forward ir angles. They attributed this 
increase over the free jtN value of 9 to a reduction in the 7 r“p 
quasifree amplitude through destructive interference with another 
process. They identified the interfering process as aN knockout. 

The present study includes features of both types of 
experiment by detecting prompt de-excitation 7 -rays in coinci- 
dence with the outgoing pions or protons. It thus combines 
kinematic information with the ability to measure excitation of 
specific nuclear states. 

We studied 200 MeV ir'*' incident on a target. Outgoing 

charged particles (scattered and knocked out protons) were 
detected by scintillator telescopes at 30*, 60*, 90*, 120 * and 
150* and identified by their dE/dX. Coincident nuclear 7 -rays 
were detected in one of two Ge(Li) detectors and identified by 
their energy. One can thus study several different (jr,jrN) 
reaction channels. For example, -y-rays in coincidence with 

a or a p can result only from a direct proton knockout with no 
charge exchange. 23jjg -y-xrays in coincidence with a can result 
only from a direct neutron knockout. On the other hand 23j^g 
7 -rays in coincidence with a proton can only result from a charge 
exchange reaction. 

24ng was chosen as a target because single nucleon removal 
from 24j^g results in mirror nuclei 23j^g ^^d 23jja. This provides 
a test of both single proton and single neutron removal 
mechanisms with pions. Furthermore, the single nucleon removal 
spectroscopic strengths ( 2 dg ^2 ^^ 1 / 2 ^ from 
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24j^g 

are concentrated in two low-lying excited states which 
7 -decay directly to the ground state [ 6 ]. For both ^^d 

^^Mg, the first excited states («0.15 MeV, 5/2'*") have spectro- 
scopic factors of -4 to 6 , and the 1/2" excited states (at 2.64 
MeV in 23jia and at 2.77 MeV in 23jjgj have spectroscopic factors 
of «4 as determined from the analysis of single-nucleon removal 
reactions on 24j^g jgj. This yields for the S/2'*' levels an 
occupation number C^S » 2,Sone-half the shell limit of 4 in 

^'^Mg. (C is the isospin coupling coefficient, (Tfrf | l/2rfilTiri) 
s yi/2.) Other bound excited states have considerably smaller 
spectroscopic factors. They predominantly feed the first excited 
state, but the combined effect of this 7 -ray feeding should be 
less than *25% of the overall strength of the first., excited state 
(if all states were populated_ in proportion to their spectro- 
scopic factors) . 

The suitability of a target for (jr,wN) reaction work 

was established previously by the results of an inclusive study 
of 7 -rays from jt'*’ reactions on ^^Mg [3]. The two above-mentioned 
states in 23ua and 23j^g with large spectroscopic factors were 
strongly excited; there was no evidence for background 7 -rays 
that might overlap these states. Another reason for the selec- 
tion of a 24j4g target was its suitability for a parallel study 
[7] of the angular correlation of 7 -rays from (ir,w') scattering. 

A feasibility study of the techniques employed in this work 
was undertaken at LAMPF using a single gamma-ray detector in 
coincidence with a single charged particle telescope [ 8 ]. The 
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results of that study suggested the need to develop a large-scale 
coincidence measurement system sensitive to de-excitation gamma 
rays, knockout nucleons and scattered pions. Such an improved 
system was developed and used in the present work. 


II. EXPERIMENTAL APPARATUS AND PROCEDURES 

The experiment was performed with a 200 MeV beam from the 
high-energy pion channel (P^) of LAMPF. This beam had a contami- 
nation of 6 . 6 % m'*’ and 2 . 0 % e"*" as well as a muon halo of roughly 
three times the beam diameter. The beam spot size was typically 
2.5 cm in diameter, arid the momentum resolution was =0.5%. The 
target consisted of natural magnesium metal (79% with an 

average density of 0.57 ± 0.02 g/cm^. 

The experimental geometry was defined by six scintillation 
telescopes for charged particle detection and two Ge(Li) spectro- 
meters to detect 7 -rays (see Fig. 1) . Each of the six particle 
telescopes consisted of the six NE 102 scintillators (Q,a,E, 
rear, left-side veto, and right-side veto) . Each scintillator 
was coupled to a 5 cm photomultiplier tube, except for the E 
scintillator, which was coupled to two 12.5 cm photomultiplier 
tubes, one at each end. The detector thicknesses were 0.160 cm 
(Q) , 0.320 cm (A), 0.635 cm (vetos and rear scintillators) and 
15.75 cm (E) . The dimensions of each telescope component were 
the same for each telescope with the exception of the Q scintil- 
lators. The n scintillators for telescopes 1, 5 and 6 at ± 

30* and 150* had to be moved further from the target to avoid the 
beam halo; they were made correspondingly larger so that all the 
telescopes subtended the same solid angle. 

The n and a counters together defined the solid angle ("0.18 
sr) of each telescope, and the E scintillator determined the 
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particle energy. All three scintillators were used for particles 
identification. The rear scintillator tagged particles which 
had not stopped in the previous scintillators. The two side veto 
scintillators tagged particles scattering out of or into sides of 
the E scintillators. During off-line data analysis, however, 
this feature was not used. Calculation of the telescope solid 
angles was performed following the method of Gotch and Yogi [9], 
considering the size and location of the n, A and E scintillators 
for each telescope. All telescopes except number 6 were mounted 
together on a pivoting table with their axis directly under the 
target centerline. Telescope 6 was mounted on a similar but 
smaller table pivoting on the same axis. 

Two Ge(Li) 7 -ray spectrometers were used in the present 
experiment, an Ortec 9% efficient and a Princeton Gamma-Tech 11% 
efficient lithium drifted germanium detector. Both detectors 
were fitted with NE 102 anti-coincidence scintillators cups to 
tag charged particles entering them. They were mounted on one 
rolling table to facilitate positioning and shielding. One 
detector was located at -75* and the other at -122* relative to 
the pion beam line axis (see Fig. 1) . Additional experimental 
details are given in Ref. [7]. 

The beam intensity was monitored with a 7.6 cm thick ion 
chamber filled with argon, and the beam profile was monitored 
with a LAMPF wire chamber system [10]. The absolute cross 
sections were normalized to the differential inelastic n scatter- 
ing cross sections from the 2‘*‘ state of [11]. 
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The six telescopes were calibrated in energy by tuning 
the channel for low-intensity protons at 50 MeV, 133 MeV and 
191 MeV and placing each of the telescopes in the beam. The 
telescopes were also calibrated with the pion and proton quasi- 
elastic scattering peaks from the experimental runs as well as 
with the maximum energy deposited in the E scintillators for 
pions and protons. The telescope calibration runs were also used 
to determine the efficiencies of the six telescopes. They were 
found to average 96 ± 2%. 

The energy responses of the two Ge(Li) detectors were 
periodically calibrated by placing ^^®Th, ®^Mn and ^^*^Cs sources 
at the target location. Well-known strong 7 -ray peaks in the 
experimental data provided additional energy calibration, 
including the first excited state to ground transition and 

the 23jja first excited state to ground transition. The maximum 
deviation of the calibration data from a linear fit was 0.7 keV. 
Relative and absolute Ge(Li) detector efficiencies were also 
determined in the source calibration runs. 

A valid data event consisted of a coincidence between a 
particle telescope signal and a Ge(Li) 7 -ray. For each event, 
pulse heights were digitized for n. A, E (two photomultipliers) , 
and rear scintillators as well as the Ge(Li) detectors. All data 
were read into a PDP-11 computer and written on magnetic tape for 
later replay. 

Particles were identified from their a and E pulse heights 
using the method of Goulding et al . [12]. A particle which 
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passed completely through the E scintillator (75 MeV pions) and 
160 MeV protons) as indicated by the rear scintillator was 
treated using the method of England [13]. Figure 2 shows a 
typical dE/dx vs. E dot plot (for Telescope 1, 30*), with the 
pions and protons identified. 

During off-line data analysis, spectra were accumulated for 
7 -rays in one of the two Ge(Li) detectors that were in coinci- 
dence with pions or protons in any one of the six particle 
telescopes. Statistics were not sufficient to allow meaningful 
cuts on pion or proton energy. Figure 3 shows the spectrum in 
the Ge(Li) detector at 75* in coincidence with a pion or a proton 
in any one of the six telescopes. Figure 4 shows the random (off 
timing peak) spectrum for this detector. 

The 7 -rays in these spectra were identified by their energy, 
and areas were determined by sTimming channels and subtracting 
background. Cross sections were calculated from the relative 
areas and from the Ge(Li) and the particle telescope efficien- 
cies. As noted above, the cross sections were normalized to the 
2^Mg 2’*' differential inelastic scattering cross sections of 
Bolger [11]. Major sources of error were the following: 
statistical errors, errors in the absolute normalization to the 
inelastic data («15%) , Ge(Li) efficiency calibration (*9%) , and 
telescope solid angle determination (»6%) . 

In addition to the strong first excited t/2'^ states at 0.439 
MeV in 23jja and 0.450 MeV in (see Fig. 3) , there was 

evidence for the fourth excited 1/2“ state in ^^Na at 2.64 MeV; 



but its Doppler-broadened width prevented a differential cross 
section measurement. There was no sign of its mirror state in 
^^Mg. 


III. COMPARISON OF ANGULAR DISTRIBUTIONS WITH REACTION MODELS 

The experimental differential cross sections for production 
of the 0.439 and 0.450 MeV 7 -rays in coincidence with outgoing 
pions or protons are listed in Table I for the two Ge(Li) 
detectors. The averages (last column of Table I) were calculated 
using the inverse of the fractional errors as weights. Where two 
cross sections differed greatly, the errors were increased in 
order to be more conservative. Data from Telescopes 1 and 6 were 
averaged to give one data point at 30*. 

The cross sections were extracted from the data by assuming 
isotropic 7 -ray correlation with the outgoing pion or proton. 

This assumption would be rigorously true for direct plan wave 
nucleon knockout and can be seen to be approximately true within 
experimental uncertainties by comparing the relative cross 
sections of Ge(Li) 1 and 2 listed in Columns 4 and 5 of Table I. 
(Ge(Li) 1 and 2 had an angular separation of 50*.) This absence 
of angular correlation is in contrast to the expected strong 
sin^a^^q correlation that was observed in this experiment for 
inelastic jt"*’ scattering to the 1.37 MeV first excited state of 
^^Mg [ 7 ]. (g^q is the angle between the 7 -ray and the inelastic 
momentum transfer direction.) 

Proton and pion differential cross sections were compared 
with predictions of an intranuclear cascade (INC) code developed 
by Frankel et al . [14] and with the predictions of a simple 
plane-wave impulse approximation (PWIA) which assumes only a 
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single collision of the incident pion with a nucleon. In Section 
IV, cross section charge ratios are compared with these models 
and also with a model that assumes final-state nucleon charge 
exchange (NCX) [15]. 

The INC predictions were based on 5*10'^ cascades of 200 MeV 
w'*’ on 24ug. The program output was sorted to yield differential 
and proton cross sections at the angles measured in the 
present experiment for events in which the final nucleus was ^^Na 
or ^^Mg in a bound state. In performing these calculations, that 
part of the code that evaluates evaporation subsequent to the 
cascade was not used. Instead, a ^3Na or ^^Mg nucleus was 
assumed to retain its identity if its excitation energy calcu- 
lated by INC following the cascade was less than its known 
particle stability energy. (The INC code does not include 
details of nuclear structure or predict specific nuclear states. 
It has, hwoever, no free parameters.) Since the INC-calculated 
contributions to the 23jja/23ng bound states arose from the entire 
nuclear volume and not just from the nuclear surface (i.e. 1 ^ 5/2 
shell nucleons), one should multiply the. INC results by the ratio 
r » SpS”(jP) •F(J^,5/2"'')/SpS"(JP) , where S“(J^) is the spectro- 
scopic factor for neutron or proton removal to a 23^a or 23j{g 
bound state of spin parity and F(J^, 5/2"^) is the relative 
fraction of 7 feeding from an initial state down to the 5/2'*' 
first excited state. Using spectroscopic factors and 7 branching 
ratios from ref. [6], r = 0.6 for both 23jja and ^^Mg. The INC 
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results quoted in the present paper have been reduced by this 
factor. 

Analysis of the output from the INC calculation showed that 
a large fraction («90%) of production of a bound a state of the 
A-*l residual nuclei results from a single collision of the 
incident pion pion on a nucleon with no further interaction. The 
probability of the residual nucleus being left with a given 
excitation decreases monotonically with increasing excitation 
energy. 

The PWIA predictions for the cross sections were calculated 
using the semiempirical free wN phase shift of Rowe et al . [16]. 
The resulting cross sections were reduced, at small pion scatter- 
ing angles, by Pauli blocking using a degenerate Fermi sphere 
uniformly filled up to a momentum of kp * 270 MeV/c. The Pauli 
blocking was calculated from the quantity (V-N)/V, where N is the 
phase space volxime common to the Fermi sphere and a similar Fermi 
sphere whose center is displaced by an amount q * 2k^ sin 9 ^/ 2 , 
and V is the sphere's volume. This caused the resulting cross 
section at « 0* to be 0 and at - 80* to approach the free 
ir-N cross section. This cross section was then multiplied by the 
ldg ^2 proton or neutron occupation number for the first excited 
state, = 2.2 defined above. It was necessary to addition- 

ally scale the resulting PWIA prediction to the data by, multi- 
plying by factors of -1/3 to 1/5. This indicates the 
predominance of other processes such as pion absorption and 
secondary scatterings. 
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The PWIA calculation was performed only to estimate the 
trends of the experimental data. It, however, demonstrates the 
role of Pauli blocking in reducing the forward angle pion 
scattering from the free ttN. 

A. Pion angular distributions 

Figure 5 shows experimental and calculated angular distribu- 
tions for outgoing pions that are in coincidence with 7 rays from 
the first excited states of 23j^g 23 solid curve 

represents the Pauli-blocked PWIA results described above; the 
open circles are the results from the INC calculation. Without 
Pauli blocking, the pion differential cross section would rise 
steadily as 8 decreases below -60*. At = 30*, the PWIA cross 
section is reduced by a factor of -2 relative to the free case. 
Both INC and PWIA calculations (and the data) display Pauli 
blocking with decreasing 8j the INC cross section falls off more 
rapidly than the PWIA cross section. An effect that could 
account for this discrepancy is nuclear shadowing of forward- 
scattered pions for the INC calculation. In the PWIA calcula- 
tion, this effect has not been included. 

B. Proton angular distributions 

Figure 6 displays the experimental angular distributions for 
outgoing protons in coincidence with 7 -rays from the first 
excited states of 23j^g 23jja^ together with PWIA and INC 
results (solid curves and open circles, respectively) . The 
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angular distributions have the general shape of free w-N scatter- 
ing, in which case no protons would be emitted at angles greater 
than 90* . 

The E scintillator thickness (15 gm/cm^) was insufficient to 
permit derivation of pion energy spectra but was adequate for 
determination of proton energy spectra by use of a fold-back 
procedure [7]. The 30* and 60* 7 -coincident proton spectra are 
shown in Fig. 7. The arrows indicate the energies for free wN 
scattering; the dashed lines indicate the INC results. Coinci- 
dent proton spectra for flp a 90* had few total counts; as 
expected, no peak was observed. The apparent differences between 
the 60* spectra and the INC predictions may be instinimental ; the 
steep left-hand shoulder is due to electronic low-energy cutoff, 
and the high-energy tail is partially due to scintillator energy 
resolution. The cross sections at 60* were corrected for the 
low-energy cutoff. 



IV. ABSOLUTE CROSS SECTIONS 


Table II compares experimental and calculated angle- 
integrated absolute cross sections. It shows that the observed 
cross sections are -1/3 to -1/5 of those predicted by the PWIA 
but are 3 times those predicted by INC. Particularly at the a 
resonance energy, the assumption of a single jtN interaction which 
is implicit in the PWIA is questionable. The 200 mb free jtN 
cross section at the resonance yields a spatial width of -5f 
resulting in a "swelling” of the nucleon to a size encompassing 
up to two adjacent nucleons. This "swelling" yields total n- 
nucleus cross sections that are typically twice the nuclear 
geometrical cross section. For example, in ^*^A1, the geometrical 
cross section is -450 mb whereas the total pion cross section is 
-960 mb [16] and is divided approximately equally among pion 
capture, nuclear elastic and non-elastic processes. The cross 
sections that remain for production of states observed in the 
present experiment are thus a small fraction of the total cross 
sections and would be sensitive to variations in any of the above 
components . 

The ratios a„(^^Na)/ap(23fig) for the coincident differen- 
tial cross sections have values of -4, approximately independent 
of 9^ (see Table 1) . This is in disagreement with the results of 
Kyle et al . [5] at T^ = 240 MeV in which the corresponding 

charge ratio R = a(jr'*’)/<7(«“) for ^®0(ir^,7r^p) reaches a 

^ • s • 

very large value, R > 30, for forward angles, 9„ < 35*. Kyle gt 
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al . suggests that this enhancement over the quasifree value of 9 
comes from a reduction in the ir“p cross section, as w'^'p enhance- 
ment is unlikely. In our experiment, however, the angle- 
integrated coincident 5/2'*’ cross section, which corre- 

sponds to the cross section of Kyle et al . , is relative to 
PWIA, the largest of all four measured cross sections (see Table 
II, column 7) . It may be noted that the enhancement found by 
Kyle et al . occurs at scattering angles where Pauli blocking is 
predominant . 

Our and p coincident 23jja cross sections are in agreement 
with equivalent cross sections obtained from <aata of 

Piasetzky et al . [4] who used a doxible arm spectrometer system. 
Their cross section values approximately equal those for free ir'^p 
scattering; and ass\aming that four p shell protons are available, 
their effective participation ratio is approximately 0.25, which 
agrees with our results for ^^Na (Table II, column 7) . However, 
their cross section is -60% greater than our equiva- 

lent coincident ^^Kg cross section. 

The INC predictions as mentioned above fall below the 
experimental results by a factor of «l/3. An examination of the 
particle histories generated by the INC code yields the following 
further information on the 24Mg(,^„p) reaction: 

A. INC predicts total, elastic, and absorption cross 
sections of 980, 400, and 230 mb, respectively, in 
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general agreement with 960, 380 and 218 mb, respec- 
tively, as measured by Ashery et al . [17] for 245 MeV 
n* + 27 ^ 1 . 

B. INC indicates that 75% of the a's decay before striking 
a nucleon. This is due to the short free decay length 
of the A (*0.4 f ) . Our PWIA calculation indicates that 
Pauli blocking of the decay is not predominant - only 
"1/3 of the A decays are Pauli-blocked. 

C. According to the INC results, pions from A decay 
predominately do not escape the nucleus; 75% of the 
pions from a decay strike a nucleon to re-form another 
A. Using the a decay probability given above in B, 
this yields an average of two sequential a's formed for 
each T » 200 MeV pion incident on the nucleus. In a 
single scattering, the pion loses an average of 60 MeV 
lada kinetic energy; after two or more pion scatterings 
through A formation, will have dropped considerably 
below resonance energy. INC yields a pion scattering 

mean free path of 1.2 f; cf. = ^/p^n total 
» 2.3 f. (Pion absorption occurs only through aN-»NN) . 

D. INC indicates that approximately 70% of the protons 
resulting from a decay escape the nucleus without 
further scattering. INC yields Xp » 4.5 f; this in 
agreement with Xp = l/p<rp > 5 f obtained from free, but 
Pauli blocked [18] NN total cross sections [19] and a 
2^Mg uniform nuclear matter density p of radius 1.315 f 
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[20]. The INC result is also in agreement with 
estiamtes by Schiffer [21]. Of the scattered a decay 
prbtons, approximately half undergo charge exchange 
before escaping in the INC calcualtions . Since only 
30% of the protons do not escape, this yields a proba- 
bility for NCX of »15%. 

E. Approximately 40% of the total cross section for 
incident pions results in pion capture (aN-*NN) , 
according to the INC calculation. Measurements by 
Ashery et al . [17] indicate a ratio of capture-to-total 
cross section of s30%. 

In a recent paper, Frankel et al . [22] explored the effects 
of using a different nucleon momentum distribution in the INC 
code. The usual version (the one used in the present paper) 
assumes a local, degenerate Fermi gas distribution. This was 
changed to a shell model with harmonic oscillator wave functions. 
The latter momentum distribution improved the INC code agreement 
with the (ir,»rp) angular correlation data of Piasetzky et al . [4]. 
Since the jr“ induced nucleon removal is likely to be a surface 
reaction, the use of a more realistic momentum distribution could 
improve agreement with the present results. 

The INC results discussed above suggest that pion multiple 
scattering, occurring mainly by sequential a production and 
decay, is a predominant process, being more important 
(x„ » 1.2 f) than nucleon multiple scattering (Ap » 5 f) which 
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was proposed some time ago [15] as a major process in pion- 
induced nucleon knockout at ^ resonance energies. In the latter 
process the nucleon from a decay undergoes subsequent incoherent 
nuclear scattering with a probability of nucleon charge exchange 
(NCX) , P = 0.1-»0.2 [15] as determined from cross section ratios 
obtained in a[^ 2 c(„+^X) llc]/a[ 12 c(ff“,X) ^^C] activation 
experiments [ 1 ] . In the present experiment , both jt'*'- and p 
coincident cross sections for de-excitation 7 -rays are 
determined. Hence the present experiment provides a more 
sensitive test of the NCX model than the previous activation 
experiments, which sum over these two final states. 

The NCX calculation predicts cross sections for the final 
states 2 ^Mg + , 23jig + 23jja + jr'*’, and ^3jja + p to be in the 

ratio of (1+9P) :2: (9+P) : (9+P) . Since the 23^3 residual nucleus 
must be accompanied by both and p, the final states 23jja + ir'*’ 
and 23jja +p must have the same cross sections regardless of the 
reaction model. Hence, the two experimental 23jia cross sections 
(see Table II) were averaged, yielding a[23Na(Av)] *» 48 ± 8 mb. 
There are then only two independent cross section ratios. Let 
them be the ratios of and cr[23Mg+p] to a[23Na(Av)]. 

The first ratio, 

Rl - <7[23Mg + 7r+]/<7[23Na(Av) ] = .23 ± .04, 

yields P = .12 ± .04, consistent with values of P obtained from 
activation work. (Quasifree =» 1/9) . However, the second 
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ratio, obtained from the proton component of the final 

state, has a value 

R 2 * cr[23Mg + p]/a[23Na (Av) ] « .33 ± .05, 

yielding a negative value of P, P = -2.9 ± .9, which is inconsis- 
tent with NCX, even though the uncertainty is quite large. 
(Quasifree R 2 = 2/9). Because of electronic cutoff and the lack 
of scintillators in the forward direction where the proton cross 
section peaks, we may be undercounting proton events and hence 
R 2 « However, increase in R 2 would worsen agreement with NCX. 

A comparison of our experimental results with the results of 
activation measurements can be made by sxunming the measured 
and p components of the ^^Hg cross section (see Table II) . This 
yields a ratio of 

a[(23Mg + 7r+) + (23Mg + p) ]/a[23Na(Av) ] » .56 ± .06, 

which in turn yields a value of P » .24 ± 0.07, consistent with 
the value of P obtained from activation measurements. 

The above results indicate that whereas the NCX model can 
explain inclusive cross section results like those obtained by 
activation measurements, it is inconsistent with the more 
exclusive cross sections obtained in the present experiment. In 
a recent paper, Ohkubo and Liu [23] include the effects of 
quantum mechanical interference between quasifree and non- 
quasifree (NCX and k charge exchange) reaction processes using 
distorted waves. Their calculations result in significantly 
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better agreement with the experimental results for 
cross section ratios [1] than the previously incoherent NCX 
calculations [15]. In asubsequent paper Ohkubo, Liu et al . [2] 
conclude that both NCX and the interference effects decrease 
considerably in magnitude as the target mass is increased from 
A » 12. Their results suggest that these effects are small for 
an A = 24 target. 

A process in which an initial a siibsequently interacts with 
a nucleon in a (AN) state can, by itself, reproduce measured 
values of and R 2 . Such a process, in which one of the (AN) 
decay nucleons subsequently remains in the nucleus, yields 
Rl = 0.22, R 2 = 0.37, i.e. values close to those observed. 
Although there has been evidence for a possible (aN) attrac- 
tive potential [24], an examination of the magnitude of pion 
double charge exchange cross sections casts doubt on this 
process. Even after allowing for the isospin recoupling, which 
yields for pion double charge exchange a (AN) component *1/3, 
any aN contribution that is sufficiently large to yield a 
reasonable (w,jrN) reaction would result in a pion double charge 
exchange cross section too high by at least a factor of 10. 


V. CONCLUSION 


The experimentally determined pion and proton differential 
cross sections for de-excitation 7 -rays from the 5/2'*’ first 
excited states of and in coincidence with outgoing, 

or p from 24Mg(^+^„N) have been measured and compared with the 
results of calculations based on several (ir^wN) reaction models, 
in particular on a Pauli blocked plane-wave impulse approximation 
(PWIA) , on intranuclear cascade (INC) models, and on charge 
exchange of the outgoing nucleon (NCX) . 

Both the PWIA and the INC calculations reproduce the 
approximate shape of the observed ir"*" and p angular distributions, 
but not the magnitude. The PWIA calculation, which was done 
primarily to illustrate the trends in the data, resulted in cross 
sections that were a factor of -3 to 5 too large. On the other 
hand, the INC calculations, which can be considered absolute, 
were a factor of *3 too small. These calculations indicate that 
rescattering of the outgoing pions is a more important process 
than interaction of the outgoing nucleons. The NCX model is put 
to a more sensitive test by the present experiment than by 
previous activation experiments, since and p coincident cross 
section ratios for de-excitation 7 -rays are determined separately 
rather than together. The NCX results are inconsistent with 
experimental results. 

These comparisons with several reaction models suggest that 
a more detailed description of the jtN interaction in a nucleus. 
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such as the A-hole model of Hirata, Lenz and Thies [25], may be 
needed for a better understanding of the processes involved in 
the (jr,)rN) reaction. 
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FIGURE CAPTIONS 

1. Experimental Geometry. The rear and side veto counters 
surrounding each E counter are shown but not labeled. 

2. A dE/dx vs. E dot plot for telescope 1 at 30*. The E signal 
was from one of the two photomultipliers of that 
scintillator. 

3. Ge(Li) #1 7 -ray spectrum in coincidence with a or p from 
24Mg(jr‘^,irN) in any one of the six particle telescopes (low- 
energy portion) . 

4. Ge(Li) #1 7 -ray spectrum in random coincidence (low-energy 
portion) . 

5. Differential cross sections of outgoing from 24jjg(^+^^+jl) 
in coincidence with 7 -rays from the first excited states of 
^^Na and 23j^g compared with Pauli-blocked plane-wave impulse 
approximation (PWIA) and intranuclear cascade (INC) calcula- 
tions. PWIA values have been multiplied by 0.20 for 23jja 
and by 0.38 for 23jjg. 

6 . Differential cross sections of protons from ^^Mg(ir'^,pp) in 
coincidence with 7 -rays from the first excited states of 
^^Na and ^3ng compared with Pauli-blocked plane-wave impulse 
(PWIA) and intranuclear cascade (INC) calculations. PWIA 
values have been multiplied by 0.17 for 23jja and by 0.28 for 
2^Mg. PWIA and INC values at backward angles are <0.1 
mb/sr and hence do not show up on the semi-log plot. 



7 . Energy spectra of protons from 24jjg(^+,p) detected at 30 ' 
and 60 * in coincidence with 7 -rays from the first excited 
states of 23jja and compared with INC calculations 

(dashed line) . The arrows indicate the energies for free 
wN scattering. 



Table I 


Experimental differential cross sections for (tt, jtN) . 


ff^(23Na) is the differential cross section for production of the 
2^Na first excited state (5/2"^/ 0.439 MeV) in coincidence with a 
Similar definitions apply for the other cross sections. 
Results are shown for each Ge(Li) and as an average which was 
weighted by the fractional errors. Telescope 1 was averaged with 
Telescope 6 for the 30* results. 


Reaction 

Telescope 

Angle 

Ge(Li)l H 

Ge(Li}2 

|| Average || 

«r,(23Na) 

1 

30* 

7.3±1.9 

4.611.2 

5.211.1 


6 

30* 

4.8±1.3 

1.911.3 



2 

60* 

4.4±1.2 

3.310.9 

3.910.9 


3 

90* 

4.011.1 

1.610.5 

3.011.0 


4 

120* 

6.311.7 

4.111.1 

5.211.0 


5 

150* ■ 

6.211.6 

4.211.1 

5.211.0 

a,(23Mg) 

1 

30* 

1.510.5 

1.910.5 

1.410.5 


6 

30* 

1.410.5 

0.3110.15 



2 

60* 

1.110.4 

0.7410.25 

0.9010.22 


3 

90* 

0.8210.34 

0.9710. 30 

0.9210.23 


4 

120’ 

1.210.5 

0.8910.29 

1.010.25 


5 

150* 

0.7210.28 

1.310.4 

1.110.30 , 

<^p(2^Na) 

1 

30* 

15.014.0 

16.014.0 

15 .014.0 


6 

30* 

18 .014.0 

11.013.0 



2 

60* 

4.811.3 

4.011.1 

4.410.9 


3 

90* 

0.5910.28 

0.8510.30 

0.7610.20 


4 

120* 

1.510.5 

0.8710.31 

1.210.3 


5 

150* 

1.310.5 

0.4510.20 

0.9310.27 

<7p(23Na) 

1 

30* 

7.612.0 

7.711.9 

5.811.4 

6 

30* 

5.111.4 

2 . 810 .80 



2 

60* 

2.110.7 

1.310.4 

1.710.4 


3 

90* 

0.4810.24 

0.3110.15 

0.3910.14 


4 

120* 

0.4510.25 

0.5910.25 

0.5410.18 


5 

150* 

0.9310.38 

0.2810.15 

0.6910.22 
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